Reverse undercompressive shock structures in driven thin film flow 
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We show experimental evidence of a new structure involving an undercompressive and reverse 
undercompressive shock for draining films driven by a surface tension gradient against gravity. The 
reverse undercompressive shock is unstable to transverse perturbations while the leading undercom- 
pressive shock is stable. Depending on the pinch-off film thickness, as controlled by the meniscus, 
either a trailing rarefaction wave or a compressive shock separates from the reverse undercompressive 
shock. 
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We consider a thin fluid film wetting a solid substrate 
and driven by a thermal gradient against gravity. This 
system exhibits a range of phenomena including two 
kinds of shocks (Lax and undercompressive) and a fin- 
gering instability that forms from a flat interface. A key 
issue that we explore is the role played by the meniscus 
at the lower boundary. 

Insight is provided from two previous studies. In exper- 
iments by Ludviksson and Lightfoot Q (LL) , a partially 
immersed substrate was pulled from a liquid bath and 
allowed to drain. A temperature gradient drove the film 
up the substrate against gravity; no fingering instabil- 
ity was reported. Cazabat et al. (Cal) created very thin 
capillary driven films |^,[| by leaving the lower end of the 
substrate in a liquid reservoir; the resulting meniscus con- 
trolled the climbing film thickness. A strong Marangoni 
stress created a capillary ridge, resulting in a fingering 
instability similar to that in gravitationally driven films 

In both types of experiments, the film thickness is a 
key parameter. In the Cal configuration, inclining the 
substrate changes the curvature of the meniscus, result- 
ing in a range of film thicknesses || . Cazabat et al. ||]7]] 
show that the bulk film thickness can increase more than 
an order of magnitude with decreasing inclination angle 
a from the horizontal. As a varies, a transition occurs 
from a stable flat front to a front that is unstable to fin- 
gering. 

This transition reflects a basic change in the dynam- 
ics of the film: very thin film fronts show a single shock 
profile that is unstable to fingering || while thicker films 
exhibit an undercompressive-compressive shock pair [|]]|] 
in which the leading undercompressive (UC) front is sta- 
ble |§,|l0|. The experiments in are, to our knowledge, 
the first examples of naturally occurring UC shocks for 
systems subject to a scalar conservation law. 

In this Letter we present experimental evidence of a 
new structure involving a pair of UC shocks. This was 
recently proposed by Munch [ |PT| to explain dynamics 
mentioned in the original LL work pi. Using a lubri- 



cation model for the meniscus region, he computes film 
dynamics for an initial condition corresponding to the LL 
experiment. In this model, the meniscus region pinches 
off a body of fluid that travels up the substrate almost 
as a solitary wave. This structure contains leading and 
trailing fronts, both undercompressive. The trailing un- 
dercompressive front has never been documented in ex- 
periments, although its existence was rigorously proved 
]l2[ . Munch calls this trailing front a 'reverse undercom- 
pressive' (RUC) shock because it involves a thicker film 
receding upward from a thinner film, which is the reverse 
of the type found in , in which a UC shock describes 
a thicker film advancing into a thinner region. 

We compare theory for the film dynamics to our ex- 
periments. Munch uses a full meniscus model to describe 
the film and meniscus. In his simulations, the meniscus 
quickly pinches off the draining film, and later acts as a 
boundary condition for the climbing film structure. In 
the early stages of the experiment, we clearly observe 
this pinch-off process. Here we model the dynamics of 
the draining film using a well-known lubrication approx- 
imation with a 'depth averaged' velocity 
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We model the meniscus as a boundary condition that we 
discuss later. In the above, 7 denotes the surface ten- 
sion, t = d-y/dx denotes the surface tension gradient, 
a the angle of inclination (from the horizontal) of the 
plane, g the gravitational constant. Also, x is along the 
direction of the flow. The coefficient of e x in Eq. (|l|) 
represents convection due to the surface tension gradi- 
ent and the component of gravity tangent to the sur- 
face. The component of gravity normal to the surface 
has a negligible effect. We couple Eq. ([!]) with mass 
conservation, ht + V • (hV) = 0. To understand the 
shock dynamics, we ignore perturbations transverse to 
e x and consider solutions h depending only on x and 
t: h t + (f(h)) x = ~(("{/3v)h 3 h xxx ) x . The flux satisfies 
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f(h) = (rh 2 /2 — sinapgh 3 /3)/ri. We rescale to dimen- 
sionless units as in 0: h — Hh, x = xl, and t = 
Tt, where H = I = (^) 1/3 - U , 3 7 T , ) 1/3 , 

and T = 2^-(|r7pgsma) 1 / 3 . Dropping the * gives the 
dimcnsionless equation 

h t + (h 2 -h 3 ) x = -(h 3 h xxx ) x . (2) 

The experiments are carried out as follows. An ox- 
idized silicon wafer is partially dipped into a reservoir 
of silicone oil (PDMS, r/=100 cSt and 7 =0.0209 N/m at 
25°C) attached to a brass plate. The wafer is pulled out 
and clamped to the plate (Fig. |l]). We create a uniform 
temperature gradient along the plate by heating at the 
bottom with a foil heater and cooling at the top with 
a circulating bath cooler. We monitor the temperature 
along the plate via thermistors. A typical gradient, which 
is constant during an experiment, is shown in the inset 
of Fig. 1. The distance between the heater and cooler 
is 20 mm, and the width of the brass plate is 152 mm 
(effectively infinite width in the transverse direction). 




FIG. 1. The schematic diagram of the experimental setup. 
A typical temperature gradient is shown in the inset. 



We measure the thickness of the film by both stan- 
dard interferometric techniques and by a novel variation 
of this technique. The primary technique uses collimated 
light from a He-Ne laser (wavelength 632.8 nm) incident 
normally on the film. We record images with a CCD 
camera and a frame grabber. Additionally, at the end 
of a run, we shine an infrared laser (830 nm, 100 mW) 
at a point on the film. The film near the laser spot is 
reduced to essentially zero thickness due to the combi- 
nation of local heating and the induced surface tension 
gradient. We then measure the thickness at other points 
on the film by counting the interference fringes from the 
reference position (One fringe = 0.226/i m) Q|. We have 
less precision determining the absolute thickness of very 
thick regions, although the relative error between nearby 
fringes is small. This is reflected in the larger absolute 
error bar shown in the top right corner of Figs. 5 and 6. 

In the early stage of the film evolution, a stationary 
pinched-off portion appears in the meniscus, while the 



contact line climbs. Later, a RUC shock moves up from 
the pinched-off portion and a broadening RW appears 
between the meniscus and the RUC shock as in Fig. 0. 
In the experiments, the leading UC shock is stable, while 
the RUC shock becomes unstable and begins to finger as 
predicted jll} . The dimensionless fingering wavenumber 
(2ttI/\) is 0.39 ± .02 which is close to the most unstable 
wavenumber of 0.35 predicted by linear theory As 
the film evolves, a flat region (/irtjc) appears just behind 
the RUC shock. 




FIG. 2. (a) Interference fringes at t = 3200 sec. for a = 85° 
and r — 0.11 Pa. The vertical and the horizontal size are 14.4 
mm and 19.2 mm respectively. Below are sections (as indi- 
cated by the dashed line in (a)) at times (b) sec, (c) 800 
sec, (d) 1600 sec, (e) 2400 sec, and (f) 3200 sec. 



The film thickness h cq just above the meniscus is de- 
termined by balancing surface tension gradient, gravity, 
and curvature ||. Thus, the inclination angle and the 
surface tension gradient control h cq . Here, we fix dj/dx 
and we vary a to control h cq . Decreasing a, leads to in- 
creasing h cq and h cq — /irtjc- If ^oq > ^ruc> the a flat 
region with a CS replaces the RW just behind the /irtjc 
region as in Fig. ||. 

The position of the leading UC shock is linear in time, 
while the RUC shock shows transient nonlinear motion, 
consistent with the theoretical model (see below). The 
leading UC shock is the same as the one observed in the 
Cal experiments and the speed is the same as in . The 
RUC shock is always faster than the leading UC shock, 
so that the size of the bump dwindles. The RUC shock 
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position as a function of time for different inclination an- 
gles all collapse when the position is rescaled by I, and 
time is rescaled by T, as shown in Fig. 0. 




* ■■ / meniscus - ' ? 

jLi .^sUL/^'j 




FIG. 3. (a) Interference fringes at 2400 sec. for a = 45° 
and r = 0.11 Pa. The dimensions are as in Fig. ^. Sections 
(dashed region) are shown at (b) sec, (c) 600 sec, (d) 1200 
sec, (e) 1800 sec, and (f) 2400 sec. 
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FIG. 4. The position of the contact line (solid symbols) 
and the position of the reverse undercompressive shock (open 
symbols) versus time. Theory as predicted from the compu- 
tation from Fig. 5 shown as solid and dashed lines. 

To compare the experimental data with the model, we 
need boundary and initial conditions. Ahead of the con- 



tact line, we choose the simplest boundary condition con- 
sistent with complete wetting: a precursor model with 
h^b>0&sx^oo H|l3|,[l4|] . As shown previously 
H0jll| , b determines both the height and speed of the 
leading UC shock. This in turn determines the height and 
speed of the trailing reverse undercompressive wave. We 
find that b — 0.005 gives the best approximation of the 
leading edge of the advancing film in these experiments. 
In the model, we take x = to be the position of the edge 
of the meniscus. For the lower boundary condition, we 
assume the meniscus has just pinched off a film of thick- 
ness h = h cq with zero third derivative h xxx (0) = 0, i.e. 
that the meniscus enforces a zero net curvature gradient 
on the bulk film region. The initial condition is the sim- 
plest approximation of a draining film, that of constant 
film thickness hf for the range < x < Xf, where Xf is 
given by the initial dipping and draining of the film. The 
parameters Xf and h cq can be measured precisely in the 
experiment. The value hf is chosen to produce a total 
volume of fluid consistent with what is observed in the 
experiment. 

We numerically integrate (||) forward in time and com- 
pare with the experimental data in Figs. 2 and 3. In both 
cases h aq and Xf are determined from experimental data 
and hf is chosen to match the volume of fluid observed 
to pinch off in the experiment. The viscosity varies with 
temperature in the experiment. We choose a viscosity for 
the time rescaling of (||) that is optimized for agreement 
with the data while consistent with the experiment. 
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FIG. 5. Comparison between PDE theory and experiment 
for data from Fig. 2a. Numerical parameters are h cq = 0.025, 
hf — 0.75, Xf = 18. The bottom figure shows a close up of 
the rarefaction wave. The bar in the top right corner is the 
absolute error for h > 0.5. The absolute error for smaller h 
is shown on the bottom figure. The relative error for h > 0.5 
equals the absolute error for the small h values. 

Fig. 5 shows the theoretical thickness profile compared 
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with the experiment from Fig. 2a, the dimensionless time 
is 140. The numerical solution evolves into a structure 
composed of a leading UC shock, trailing RUC shock, 
and expanding RW, shown closeup in the bottom por- 
tion of Fig. 5. The leading UC shock speed and height 
are determined by the precursor thickness as in previous 
studies HJ^] . The RUC shock speed and height are deter- 
mined by the height of the leading UC shock, and hence 
also determined by the precursor. Fig. 5 shows a well 
developed structure, however at early times the RUC ex- 
periences a period of adjustment before it settles into a 
traveling wave shape, while the leading UC shock forms 
almost immediately. The trailing RW expands between 
the meniscus region (set by the boundary condition h eq 
in the model) and the RUC shock. The trailing side of 
the RUC shock has a thickness /iruc • If ^ruc is greater 
than h eq , a RW forms between the RUC shock and the 
meniscus. Otherwise a CS forms to connect the menis- 
cus pinch-off thickness h eq to the thickness /iruc- This 
CS is seen to slowly separate from the RUC shock. Fig. 
6, which compares the PDE solution with the data from 
Fig. 3a, in which h cq > /iruc- A CS forms to the left 
of the RUC shock, as shown close up in the bottom of 
Fig. 6. At this a of 45°, we are close to this transition 
point, and the experimental results are quite sensitive to 
changes in the initial draining film profile. 
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FIG. 6. Comparison between PDE theory and experiment 
for data from Fig. 3a, numerical parameters are h ca — 0.081, 
hf — 0.75, Xf — 15.4. The error bars are as in Fig. 

The experiments presented show clear evidence of 
the new UC-RUC double shock structure predicted by 



Munch. Our experimental results are compared with 
a lubrication model in which the pinch-off dynamics of 
the meniscus are incorporated into a boundary condition. 
There is good agreement between theory and experiment. 
The dynamics between the RUC shock and the meniscus 
can result in either a RW or a CS, depending on the 
pinch-off thickness compared with the RUC film thick- 
ness /iruc- Both situations are seen in the experiment 
and compared with theory. Future studies should address 
the nonlinear dynamics of the trailing RUC shock. Also 
the effect of temperature on viscosity is not considered in 
the model here but clearly plays a role in the experiment. 
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